Single-nucleus RNA Sequencing Research Articles

Central neurocytoma exhibits radial glial cell signatures with FGFR3 hypomethylation and overexpression.

We explored the genomic events underlying central neurocytoma (CN), a rare neoplasm of the central nervous system, via multiomics approaches, including whole-exome sequencing, bulk and single-nuclei RNA sequencing, and methylation sequencing. We identified FGFR3 hypomethylation leading to FGFR3 overexpression as a major event in the ontogeny of CN that affects crucial downstream events, such as aberrant PI3K-AKT activity and neuronal development pathways. Furthermore, we found similarities between CN and radial glial cells based on analyses of gene markers and CN tumor cells and postulate that CN tumorigenesis is due to dysregulation of radial glial cell differentiation into neurons. Our data demonstrate the potential role of FGFR3 as one of the leading drivers of tumorigenesis in CN.

Rescue of cardiac dysfunction during chemotherapy in acute myeloid leukaemia by blocking IL-1α.

Patients with acute myeloid leukaemia (AML) suffer from severe myocardial injury during daunorubicin (DNR)-based chemotherapy and are at high risk of cardiac mortality. The crosstalk between tumour cells and cardiomyocytes might play an important role in chemotherapy-related cardiotoxicity, but this has yet to be demonstrated. This study aimed to identify its underlying mechanism and explore potential therapeutic targets. Cardiac tissues were harvested from an AML patient after DNR-based chemotherapy and were subjected to single-nucleus RNA sequencing. Cardiac metabolism and function were evaluated in AML mice after DNR treatment by using positron emission tomography, magnetic resonance imaging, and stable-isotope tracing metabolomics. Plasma cytokines were screened in AML mice after DNR treatment. Genetically modified mice and cell lines were used to validate the central role of the identified cytokine and explore its downstream effectors. In the AML patient, disruption of cardiac metabolic homeostasis was associated with heart dysfunction after DNR-based chemotherapy. In AML mice, cardiac fatty acid utilization was attenuated, resulting in cardiac dysfunction after DNR treatment, but these phenotypes were not observed in similarly treated tumour-free mice. Furthermore, tumour cell-derived interleukin (IL)-1α was identified as a primary factor leading to DNR-induced cardiac dysfunction and administration of an anti-IL-1α neutralizing antibody could improve cardiac functions in AML mice after DNR treatment. This study revealed that crosstalk between tumour cells and cardiomyocytes during chemotherapy could disturb cardiac energy metabolism and impair heart function. IL-1α neutralizing antibody treatment is a promising strategy for alleviating chemotherapy-induced cardiotoxicity in AML patients.

Chronic intermittent hypoxia elicits distinct transcriptomic responses among neurons and oligodendrocytes within the brainstem of mice.

Chronic intermittent hypoxia (CIH) is a prevalent condition characterized by recurrent episodes of oxygen deprivation, linked to respiratory and neurological disorders. Prolonged CIH is known to have adverse effects, including endothelial dysfunction, chronic inflammation, oxidative stress, and impaired neuronal function. These factors can contribute to serious comorbidities, including metabolic disorders and cardiovascular diseases. To investigate the molecular impact of CIH, we examined male C57BL/6J mice exposed to CIH for 21 days, comparing to normoxic controls. We employed single-nucleus RNA sequencing to comprehensively examine the transcriptomic impact of CIH on key cell classes within the brainstem, specifically excitatory neurons, inhibitory neurons, and oligodendrocytes. These cell classes regulate essential physiological functions, including autonomic tone, cardiovascular control, and respiration. Through analysis of 10,995 nuclei isolated from pontine-medullary tissue, we identified seven major cell classes, further subdivided into 24 clusters. Our findings among these cell classes, revealed significant differential gene expression, underscoring their distinct responses to CIH. Notably, neurons exhibited transcriptional dysregulation of genes associated with synaptic transmission, and structural remodeling. In addition, we find dysregulated genes encoding ion channels, and inflammatory response. Concurrently, oligodendrocytes exhibited dysregulated genes associated with oxidative phosphorylation and oxidative stress. Utilizing CellChat network analysis, we uncovered CIH-dependent altered patterns of diffusible inter-cellular signaling. These insights offer a comprehensive transcriptomic cellular atlas of the pons-medulla and provide a fundamental resource for the analysis of molecular adaptations triggered by CIH.

Transcriptomic atlas of midbrain dopamine neurons uncovers differential vulnerability in a Parkinsonism lesion model

Midbrain dopamine (mDA) neurons comprise diverse cells with unique innervation targets and functions. This is illustrated by the selective sensitivity of mDA neurons of the substantia nigra compacta (SNc) in patients with Parkinson’s disease, while those in the ventral tegmental area (VTA) are relatively spared. Here, we used single nuclei RNA sequencing (snRNA-seq) of approximately 70,000 mouse midbrain cells to build a high-resolution atlas of mouse mDA neuron diversity at the molecular level. The results showed that differences between mDA neuron groups could best be understood as a continuum without sharp differences between subtypes. Thus, we assigned mDA neurons to several ‘territories’ and ‘neighborhoods’ within a shifting gene expression landscape where boundaries are gradual rather than discrete. Based on the enriched gene expression patterns of these territories and neighborhoods, we were able to localize them in the adult mouse midbrain. Moreover, because the underlying mechanisms for the variable sensitivities of diverse mDA neurons to pathological insults are not well understood, we analyzed surviving neurons after partial 6-hydroxydopamine (6-OHDA) lesions to unravel gene expression patterns that correlate with mDA neuron vulnerability and resilience. Together, this atlas provides a basis for further studies on the neurophysiological role of mDA neurons in health and disease.

Transcriptomic atlas of midbrain dopamine neurons uncovers differential vulnerability in a Parkinsonism lesion model.

Midbrain dopamine (mDA) neurons comprise diverse cells with unique innervation targets and functions. This is illustrated by the selective sensitivity of mDA neurons of the substantia nigra compacta (SNc) in patients with Parkinson's disease, while those in the ventral tegmental area (VTA) are relatively spared. Here, we used single nuclei RNA sequencing (snRNA-seq) of approximately 70,000 mouse midbrain cells to build a high-resolution atlas of mouse mDA neuron diversity at the molecular level. The results showed that differences between mDA neuron groups could best be understood as a continuum without sharp differences between subtypes. Thus, we assigned mDA neurons to several 'territories' and 'neighborhoods' within a shifting gene expression landscape where boundaries are gradual rather than discrete. Based on the enriched gene expression patterns of these territories and neighborhoods, we were able to localize them in the adult mouse midbrain. Moreover, because the underlying mechanisms for the variable sensitivities of diverse mDA neurons to pathological insults are not well understood, we analyzed surviving neurons after partial 6-hydroxydopamine (6-OHDA) lesions to unravel gene expression patterns that correlate with mDA neuron vulnerability and resilience. Together, this atlas provides a basis for further studies on the neurophysiological role of mDA neurons in health and disease.

Empagliflozin and liraglutide ameliorate HFpEF in mice via augmenting the Erbb4 signaling pathway.

Heart failure with preserved ejection fraction (HFpEF) is closely associated with metabolic derangement. Sodium glucose cotransporter-2 inhibitors (SGLT2i) and glucagon-like peptide-1 receptor agonists (GLP-1RA) exert anti-HFpEF effects, but the underlying mechanisms remain unclear. In this study, we explored the anti-HFpEF effects of empagliflozin and liraglutide and the underlying molecular mechanisms in a mouse model of HFpEF. This model was established by high-fat diet (HFD) feedingplus Nω-nitro-L-arginine methyl ester (L-NAME) treatment. The mice were treated with empagliflozin (20 mg·kg-1·d-1, i.g.) or liraglutide (0.3 mg·kg-1·d-1, i.p.) or their combination for 4 weeks. At the end of the experimental protocol, cardiac function was measured using ultrasound, then mice were euthanized and heart, liver, and kidney tissues were collected. Nuclei were isolated from frozen mouse ventricular tissue for single-nucleus RNA-sequencing (snRNA-seq). We showed that administration of empagliflozin or liraglutide alone or in combination significantly improved diastolic function, ameliorated cardiomyocyte hypertrophy and cardiac fibrosis, as well as exercise tolerance but no synergism was observed in the combination group. Furthermore, empagliflozin and/or liraglutide lowered body weight, improved glucose metabolism, lowered blood pressure, and improved liver and kidney function. After the withdrawal of empagliflozin or liraglutide for 1 week, these beneficial effects tended to diminish. The snRNA-seq analysis revealed a subcluster of myocytes, in which Erbb4 expression was down-regulated under HFpEF conditions, and restored by empagliflozin or liraglutide. Pseudo-time trajectory analysis and cell-to-cell communication studies confirmed that the Erbb4 pathway was a prominent pathway essential for both drug actions. In the HFpEF mouse model, both empagliflozin and liraglutide reversed Erbb4 down-regulation. In rat h9c2 cells, we showed that palmitic acid- or high glucose-induced changes in PKCα and/or ERK1/2 phosphorylation at least in part through Erbb4. Collectively, the single-cell atlas reveals the anti-HFpEF mechanism of empagliflozin and liraglutide, suggesting that Erbb4 pathway represents a new therapeutic target for HFpEF. Effects and mechanisms of action of empagliflozin and liraglutide in HFpEF mice. HFpEF was induced with a high-fat diet and L-NAME for 15 weeks, and treatment with empagliflozin and liraglutide improved the HFpEF phenotype. Single nucleus RNA sequencing (snRNA-seq) was used to reveal the underlying mechanism of action of empagliflozin and liraglutide.

Haploinsufficiency of phosphodiesterase 10A activates PI3K/AKT signaling independent of PTEN to induce an aggressive glioma phenotype.

Glioblastoma is universally fatal and characterized by frequent chromosomal copy number alterations harboring oncogenes and tumor suppressors. In this study, we analyzed exome-wide human glioblastoma copy number data and found that cytoband 6q27 is an independent poor prognostic marker in multiple data sets. We then combined CRISPR-Cas9 data, human spatial transcriptomic data, and human and mouse RNA sequencing data to nominate PDE10A as a potential haploinsufficient tumor suppressor in the 6q27 region. Mouse glioblastoma modeling using the RCAS/tv-a system confirmed that Pde10a suppression induced an aggressive glioma phenotype in vivo and resistance to temozolomide and radiation therapy in vitro. Cell culture analysis showed that decreased Pde10a expression led to increased PI3K/AKT signaling in a Pten-independent manner, a response blocked by selective PI3K inhibitors. Single-nucleus RNA sequencing from our mouse gliomas in vivo, in combination with cell culture validation, further showed that Pde10a suppression was associated with a proneural-to-mesenchymal transition that exhibited increased cell adhesion and decreased cell migration. Our results indicate that glioblastoma patients harboring PDE10A loss have worse outcomes and potentially increased sensitivity to PI3K inhibition.

Spatial co-transcriptomics reveals discrete stages of the arbuscular mycorrhizal symbiosis.

The symbiotic interaction of plants with arbuscular mycorrhizal (AM) fungi is ancient and widespread. Plants provide AM fungi with carbon in exchange for nutrients and water, making this interaction a prime target for crop improvement. However, plant-fungal interactions are restricted to a small subset of root cells, precluding the application of most conventional functional genomic techniques to study the molecular bases of these interactions. Here we used single-nucleus and spatial RNA sequencing to explore both Medicago truncatula and Rhizophagus irregularis transcriptomes in AM symbiosis at cellular and spatial resolution. Integrated, spatially registered single-cell maps revealed infected and uninfected plant root cell types. We observed that cortex cells exhibit distinct transcriptome profiles during different stages of colonization by AM fungi, indicating dynamic interplay between both organisms during establishment of the cellular interface enabling successful symbiosis. Our study provides insight into a symbiotic relationship of major agricultural and environmental importance and demonstrates a paradigm combining single-cell and spatial transcriptomics for the analysis of complex organismal interactions.

Cell-cycle-specific autoencoding improves cluster analysis of cycling cardiomyocytes.

Our previous analyses of cardiomyocyte single-nucleus RNA sequencing (snRNAseq) data from the hearts of fetal pigs and pigs that underwent apical resection surgery on postnatal day (P) 1 (ARP1), myocardial infarction (MI) surgery on P28 (MIP28), both ARP1 and MIP28 (ARP1MIP28), or controls (no surgical procedure or CTL) identified 10 cardiomyocyte subpopulations (clusters), one of which appeared to be primed to proliferate in response to MI. However, the clusters composed of primarily proliferating cardiomyocytes still contained non-cycling cells, and we were unable to distinguish between cardiomyocytes in different phases of the cell cycle. Here, we improved the precision of our assessments by conducting similar analyses with snRNAseq data for only the 1646 genes included under the Gene Ontology term "cell cycle." Two cardiac snRNAseq datasets, one from mice (GEO dataset number GSE130699) and one from pigs (GEO dataset number GSE185289), were evaluated via our cell-cycle-specific analytical pipeline. Cycling cells were identified via the co-expression of five proliferation markers (AURKB, MKI67, INCENP, CDCA8, and BIRC5). The cell-cycle-specific autoencoder (CSA) algorithm identified seven cardiomyocyte clusters in mouse hearts (mCM1 and mCM3-mCM8), including one prominent cluster of cycling cardiomyocytes in animals that underwent MI or Sham surgery on P1. Five cardiomyocyte clusters (pCM1, pCM3-pCM6) were identified in pig hearts, two of which (pCM1 and pCM4) displayed evidence of cell cycle activity; pCM4 was found primarily in hearts from fetal pigs, while pCM1 comprised a small proportion of cardiomyocytes in both fetal hearts and hearts from ARP1MIP28 pigs during the two weeks after MI induction, but was nearly undetectable in all other experimental groups and at all other time points. Furthermore, pseudotime trajectory analysis of snRNAseq data from fetal pig cardiomyocytes identified a pathway that began at pCM3, passed through pCM2, and ended at pCM1, and whereas pCM3 was enriched for the expression of a cell-cycle activator that regulates the G1/S phase transition (Cyclin D2), pCM2 was enriched for an S-phase regulator (CCNE2), and pCM1 was enriched for the expression of a gene that regulates the G2M phase transition and mitosis (Cyclin B2). We also identified four transcription factors (E2F8, FOXM1, GLI3, and RAD51) that were more abundantly expressed in cardiomyocytes from regenerative mouse hearts than from nonregenerative mouse hearts, from the hearts of fetal pigs than from CTL pig hearts, and from ARP1MIP28 pig hearts than from MIP28 pig hearts during the two weeks after MI induction. The CSA algorithm improved the precision of our assessments of cell cycle activity in cardiomyocyte subpopulations and enabled us to identify a trajectory across three clusters that appeared to track the onset and progression of cell-cycle activity in cardiomyocytes from fetal pigs.

NFĸB signaling drives myocardial injury via CCR2+ macrophages in a preclinical model of arrhythmogenic cardiomyopathy.

Nuclear factor kappa-B (NFκB) is activated in arrhythmogenic cardiomyopathy (ACM) patient-derived iPSC-cardiac myocytes under basal conditions and inhibition of NFκB signaling prevents disease in Dsg2mut/mut mice, a robust mouse model of ACM. Here, we used genetic approaches and single cell RNA sequencing to define the contributions of immune signaling in cardiac myocytes and macrophages in the natural progression of ACM using Dsg2mut/mut mice. We found that NFκB signaling in cardiac myocytes drives myocardial injury, contractile dysfunction, and arrhythmias in Dsg2mut/mut mice. NFκB signaling in cardiac myocytes mobilizes macrophages expressing C-C motif chemokine receptor-2 (CCR2+ cells) to affected areas within the heart, where they mediate myocardial injury and arrhythmias. Contractile dysfunction in Dsg2mut/mut mice is caused both by loss of heart muscle and negative inotropic effects of inflammation in viable muscle. Single nucleus RNA sequencing and cellular indexing of transcriptomes and epitomes (CITE-seq) studies revealed marked pro-inflammatory changes in gene expression and the cellular landscape in hearts of Dsg2mut/mut mice involving cardiac myocytes, fibroblasts and CCR2+ macrophages. Changes in gene expression in cardiac myocytes and fibroblasts in Dsg2mut/mut mice were dependent on CCR2+ macrophage recruitment to the heart. These results highlight complex mechanisms of immune injury and regulatory crosstalk between cardiac myocytes, inflammatory cells and fibroblasts in the pathogenesis of ACM.

Age-Dependent RGS5 Loss in Pericytes Induces Cardiac Dysfunction and Fibrosis.

Pericytes are capillary-associated mural cells involved in the maintenance and stability of the vascular network. Although aging is one of the main risk factors for cardiovascular disease, the consequences of aging on cardiac pericytes are unknown. In this study, we have combined single-nucleus RNA sequencing and histological analysis to determine the effects of aging on cardiac pericytes. Furthermore, we have conducted in vivo and in vitro analysis of RGS5 (regulator of G-protein signaling 5) loss of function and finally have performed pericytes-fibroblasts coculture studies to understand the effect of RGS5 deletion in pericytes on the neighboring fibroblasts. Aging reduced the pericyte area and capillary coverage in the murine heart. Single-nucleus RNA sequencing analysis further revealed that the expression of Rgs5 was reduced in cardiac pericytes from aged mice. In vivo and in vitro studies showed that the deletion of RGS5 impaired cardiac function, induced fibrosis, and morphological changes in pericytes characterized by a profibrotic gene expression signature and the expression of different ECM (extracellular matrix) components and growth factors, for example, TGFB2 and PDGFB. Indeed, culturing fibroblasts with the supernatant of RGS5-deficient pericytes induced their activation as evidenced by the increased expression of αSMA (alpha smooth muscle actin) in a TGFβ (transforming growth factor beta)2-dependent mechanism. Our results have identified RGS5 as a crucial regulator of pericyte function during cardiac aging. The deletion of RGS5 causes cardiac dysfunction and induces myocardial fibrosis, one of the hallmarks of cardiac aging.

Involution of brown adipose tissue through a Syntaxin 4 dependent pyroptosis pathway

Aging, chronic high-fat diet feeding, or housing at thermoneutrality induces brown adipose tissue (BAT) involution, a process characterized by reduction of BAT mass and function with increased lipid droplet size. Single nuclei RNA sequencing of aged mice identifies a specific brown adipocyte population of Ucp1-low cells that are pyroptotic and display a reduction in the longevity gene syntaxin 4 (Stx4a). Similar to aged brown adipocytes, Ucp1-STX4KO mice display loss of brown adipose tissue mass and thermogenic dysfunction concomitant with increased pyroptosis. Restoration of STX4 expression or suppression of pyroptosis activation protects against the decline in both mass and thermogenic activity in the aged and Ucp1-STX4KO mice. Mechanistically, STX4 deficiency reduces oxidative phosphorylation, glucose uptake, and glycolysis leading to reduced ATP levels, a known triggering signal for pyroptosis. Together, these data demonstrate an understanding of rapid brown adipocyte involution and that physiologic aging and thermogenic dysfunction result from pyroptotic signaling activation.

Co-expression of immune checkpoints in glioblastoma revealed by single-nucleus RNA sequencing and spatial transcriptomics

Glioblastoma (GBM) is one of the most malignant tumors of the central nervous system. The pattern of immune checkpoint expression in GBM remains largely unknown. We performed snRNA-Seq and spatial transcriptomic (ST) analyses on untreated GBM samples. 8 major cell types were found in both tumor and adjacent normal tissues, with variations in infiltration grade. Neoplastic cells_6 was identified in malignant cells with high expression of invasion and proliferator-related genes, and analyzed its interactions with microglia, MDM cells and T cells. Significant alterations in ligand-receptor interactions were observed, particularly between Neoplastic cells_6 and microglia, and found prominent expression of VISTA/VSIG3, suggesting a potential mechanism for evading immune system attacks. High expression of TIM-3, VISTA, PSGL-1 and VSIG-3 with similar expression patterns in GBM, may have potential as therapeutic targets. The prognostic value of VISTA expression was cross-validated in 180 glioma patients, and it was observed that patients with high VISTA expression had a poorer prognosis. In addition, multimodal cross analysis integrated SnRNA-seq and ST, revealing complex intracellular communication and mapping the GBM tumor microenvironment. This study reveals novel molecular characteristics of GBM, co-expression of immune checkpoints, and potential therapeutic targets, contributing to improving the understanding and treatment of GBM.

Integration of Single-nucleus and Exosome RNA Sequencing Dissected Inter-Cellular Communication and Biomarkers in Pancreatic Ductal Adenocarcinoma

BackgroundMounting evidence underscores the importance of cell communication within the tumor microenvironment, which is pivotal in tumor proliferation, invasion, and metastasis. Exosomes play a crucial role in cell-to-cell communication. Although single-cell RNA sequencing (scRNA-seq) provides insights into individual cell transcriptional characteristics, it falls short of comprehensively capturing exosome-mediated intercellular communication. MethodWe analyzed Pancreatic Ductal Adenocarcinoma (PDAC) tissues, separating supernatant and precipitate for exosome purification and single-cell nucleus suspension. We then constructed Single-nucleus RNA sequencing (snRNA-seq) and small RNA-seq libraries from these components. Our bioinformatic analysis integrated these sequences with ligand-receptor analysis and public miRNA data to map the cell communication network. ResultsWe established intercellular communication networks using bioinformatic analysis to track exosome miRNA effects and ligand-receptor pairs. Significantly, hsa-miR-1293 emerged as a prognostic biomarker for pancreatic cancer, linked to immune evasion, increased myeloid-derived suppressor cells, and poorer prognosis. Targeting this miRNA may enhance anti-tumor immunity and improve outcomes. ConclusionOur study offers a novel approach to constructing intercellular communication networks using snRNA-seq and exosome-small RNA sequencing. By integrating miRNA tracing with ligand-receptor analysis, we illuminate the complex interactions in the pancreatic cancer microenvironment, highlighting the pivotal role of miRNAs and identifying potential biomarkers and therapeutic targets.

Astroglial activation is exacerbated in a Down syndrome mouse model

Down syndrome (DS), also known as trisomy 21, is one of the most common chromosomal disorders associated with intellectual disability. Mouse models are valuable for mechanistic and therapeutic intervention studies. The purpose of this study was to investigate astroglial anomalies in Dp16, a widely used DS mouse model. Brain sections were prepared from one-month-old Dp16 mice and their littermates, immunostained with an anti-GFAP or anti-S100B antibody, and imaged to reconstruct astroglial morphology in three dimensions. No significant difference in the number of astrocytes was found in either the hippocampal CA1 region or cortex between Dp16 and WT mice. However, the average astroglial volume in Dp16 was significantly (P<0.05) greater than that in WT, suggesting the astroglial activation. Reanalysis of the single-nucleus RNA sequencing data indicated that the genes differentially expressed between WT and Dp16 astrocytes were associated with synapse organization and neuronal projection. In contrast, in vitro cultured neonatal astrocytes did not exhibit significant morphological changes. The expression of Gfap in in vitro cultured Dp16 astrocytes was not increased as it was in in vivo hippocampal tissue. However, after treatment with lipopolysaccharides, the inflammatory response gene IFNβ increased significantly more in Dp16 astrocytes than in WT astrocytes. Overall, our results showed that the increase in astrogliogenesis in DS was not apparent in the early life of Dp16 mice, while astrocyte activation, which may be partly caused by increased responses to inflammatory stimulation, was significant. The inflammatory response of astrocytes might be a potential therapeutic target for DS intellectual disability.

Polyploidisation pleiotropically buffers ageing in hepatocytes

Polyploidy in hepatocytes has been proposed as a genetic mechanism to buffer against transcriptional dysregulation. Here, we aim to demonstrate the role of polyploidy in modulating gene regulatory networks in hepatocytes during ageing. We performed single-nucleus RNA-sequencing in hepatocyte nuclei of different ploidy levels isolated from young and old wild-type mice. Changes in the gene expression and regulatory network were compared to three independent haploinsufficient strains for HNF4A, CEBPA or CTCF, representing non-deleterious perturbations. Phenotypic characteristics of the liver section were additionally evaluated histologically, whereas the genomic allele composition of hepatocytes was analysed by BaseScope. We observed that ageing in wild-type mice results in nuclei polyploidy and marked increase in steatosis. Haploinsufficiency of liver-specific master regulators (HFN4A or CEBPA) results in the enrichment of hepatocytes with tetraploid nuclei at a young age, affecting the genomic regulatory network, and dramatically suppressing ageing-related steatosis tissue-wide. Notably, these phenotypes are not the result of subtle disruption to liver-specific transcriptional networks, since haploinsufficiency in CTCF insulator protein resulted in the same phenotype. Further quantification of genotypes of tetraploid hepatocytes in young and old HFN4A haploinsufficient mice revealed that during ageing, tetraploid hepatocytes lead to the selection of wild-type alleles, restoring non-deleterious genetic perturbation. ConclusionsOur results suggest a model whereby polyploidisation leads to fundamentally different cell states. Polyploid conversion enables pleiotropic buffering against age-related decline via non-random allelic segregation to restore a wild-type genome.

Single cell transcriptome profiling of infrapatellar fat pad highlights the role of interstitial inflammatory fibroblasts in osteoarthritis

ObjectivesOsteoarthritis (OA) is a whole-joint disease in which the role of the infrapatellar fat pad (IFP) in its pathogenesis is unclear. Our study explored the cellular heterogeneity of IFP to understand OA and identify therapeutic targets. MethodsSingle-cell and single-nuclei RNA sequencing were used to analyze 10 IFP samples, comprising 5 from OA patients and 5 from healthy controls. Analyses included differential gene expression, enrichment, pseudotime trajectory, and cellular communication, along with comparative studies with visceral and subcutaneous fats. Key subcluster and pathways were validated using multiplex immunohistochemistry. ResultsThe scRNA-seq performed on the IFPs of the OA and control group profiled the gene expressions of over 49,674 cells belonging to 11 major cell types. We discovered that adipose stem and progenitor cells (ASPCs), contributing to the formation of both adipocytes and synovial-lining fibroblasts (SLF). Interstitial inflammatory fibroblasts (iiFBs) were a subcluster of ASPCs that exhibit notable pro-inflammatory and proliferative characteristics. We identified four adipocyte subtypes, with one subtype showing a reduced lipid synthesis ability. Furthermore, iiFBs modulated the activities of macrophages and T cells in the IFP. Compared to subcutaneous and visceral adipose tissues, iiFBs represented a distinctive subpopulation of ASPCs in IFP that regulated cartilage proliferation through the MK pathway. ConclusionThis study presents a comprehensive single-cell transcriptomic atlas of IFP, uncovering its complex cellular landscape and potential impact on OA progression. Our findings highlight the role of iiFBs in OA, especially through MK pathway, opening new avenues for understanding OA pathogenesis and developing novel targeted therapies.

Abstract 1076: Single nucleus RNA-sequencing using normal human salivary glands revealed a potential etiologic relationship between each salivary gland cell and salivary gland malignancies

Abstract Salivary glands are composed of several types of cells, and each cell type is predicted to be involved in the carcinogenesis of different types of cancers including Adenoid cystic carcinoma (ACC), Acinic cell carcinoma (AciCC), salivary duct carcinoma (SDC), myoepithelial carcinoma (MECA) and other histology. In this study, we performed single nucleus RNA-seq on three human salivary gland samples to clarify the gene expression profile of each complex cellular component of the salivary glands and related these expression patterns to expression found in salivary gland cancers (SGC) to infer cell of origin. By single nucleus RNA-seq, total 13,643 salivary gland cells were stratified into four clusters: acinar cells, ductal cells 1, ductal cells 2, and myoepithelial cells/stromal cells, with differentially expressed genes of each group. The localization of each cell group was anatomically validated by IHC of each cluster marker gene, and one group of ductal cells was found to represent intercalated ductal cells labeled with HES1. Gene ontology analysis for this group included “epithelial cell differentiation” confirming the intercalated duct. Furthermore, to elucidate the potential etiologic relationship between normal salivary glands cell type and salivary gland carcinoma histology, we compared the expression of each cluster marker to SGC public RNA-seq data. In comparison between AciCC and ACC, we used normalized RNA-seq data of 54 ACC and 21 AciCC. ACC had significantly less expression of acinar cells markers (STATH, CST2, and LPO) compared to other markers expression(***P = 4 × 10−51), unlike AciCC (P = 0.54), suggesting a possible contribution of acinar cells to the carcinogenesis of AciCC. In another cohort, we compared 16 SDC, and 38 MECA with publicly available RNA-seq data. Remarkably, HES1 and ELF3, each of two different ductal cell markers, were significantly upregulated in SDC consistent with origin of SDC from ductal cells (P=0.001, P=0.0007, respectively). Consistent with bulk RNA-seq results, expression of HES1 in IHC using clinical human SDC and MECA samples was higher in SDC than MECA, suggesting that the expression of the derived cells is involved in salivary gland carcinogenesis. Cell type expressions in specific SGC histology are similar to those found in normal salivary gland populations, indicating a potential etiologic relationship. On the other hand, there was no significant difference in myoepithelial marker ACTA2 expression between SDC and MECA, suggesting that the markers of origin are not always highly expressed, depending on cancer subtypes. Although publicly available SGC RNA-seq data used in this study are limited, the results remarkably reflect the diverse modes of carcinogenesis of SGC. These findings might provide a better understanding of the molecular basis of SGC and treatment strategies for SGC. Citation Format: Takuya Nakagawa, Chanond Nasamran, Jingjing Hu, Prakriti Sen, Theresa Guo, Kathleen Fisch, Joseph Califano. Single nucleus RNA-sequencing using normal human salivary glands revealed a potential etiologic relationship between each salivary gland cell and salivary gland malignancies [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 1076.

Abstract 6938: Identification of distinct tumor-TME ecomodules in glioma from neurofibromatosis type 1

Abstract Neurofibromatosis type 1 (NF-1) is the most common cancer predisposition syndrome in which 15-20% of affected individuals develop glioma. Large scale DNA and RNA bulk profiling showed the molecular complexity of NF-1 glioma with the tumor cellular ecosystem constituted by multiple malignant phenotypes and heterogenous immune microenvironment. However, the composition and function of infiltrating cells was hidden in the bulk tumor, and the extended granularity of NF-1 glioma tumor microenvironment (TME) remained still unexplored. Here, we collected glioma samples from 46 NF-1 patients including 22 high-grade (HGG) and 24 low-grade (LGG) tumors, and we analyzed their gene expression by single nuclei RNA sequencing. A total of 239,044 single cells were classified into tumor and non-tumor components by integrating multiple computational approaches (including genomic copy number inference, gene signature enrichment, and clustering). We defined the pattern of intra-tumor heterogeneity of NF-1 glioma cells using non-negative matrix factorization and derived 7 malignant meta-programs (MPs) that we respectively defined as Neuronal-like, EMT, Astrocyte-like, Dividing Radial Glia-like, Ependymal-like, Immune, and Glycolytic/Hypoxic-like. These MPs recapitulated normal brain cell subtypes, thereby reflecting broad cell plasticity. The non-tumor cell compartment (121,364 cells, 51%) was dissected for the characterization of the cell types that populate the TME of NF-1 glioma. We identified different subpopulations exhibiting specific immune functions within myeloid and lymphoid components. Different glioma ecomodules were highlighted by comparing the relative composition of the TME across the tumors. Recruitment and activation of cytotoxic CD8+ T cells and natural killers by an active crosstalk with dendritic and pro-inflammatory myeloid cells defined an immune-supportive phenotype that could mediate a potential anti-tumor response in low-grade NF-1 glioma (LGG immune high). Conversely, regulatory T cell infiltration and effector T cell exhaustion induced immune suppression in a low-grade glioma immune dysfunctional ecomodule. The absence of lymphocytes characterized a large set of cold tumors, mostly including high-grade glioma. Together, the complex interplay of tumor cell states with different TME compartments elucidated the existence of separate ecomodules in NF-1 glioma, with the LGG immune high TME associated with Neuronal-like and the LGG immune dysfunctional with Ependymal-like tumor cells. The Ependymal-like state also exhibits maximal association with brain-specific normal cells, including oligodendrocytes, neurons and astrocytes, whereas the HGG are enriched with Dividing Radial Glia- and Glycolytic/Hypoxic-like tumor cell states. The elucidation of different ecomodules provides novel insights for the application of targeted therapies in NF-1 glioma patients. Citation Format: Luciano Garofano, Fulvio D'Angelo, Michael Oh, Michele Ceccarelli, Franck Bielle, Marc Sanson, Anna Lasorella, Antonio Iavarone. Identification of distinct tumor-TME ecomodules in glioma from neurofibromatosis type 1 [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 6938.

Abstract 5474: Biorepository yields potential for personalized medicine approach in pediatric brain cancer patients

Abstract Childhood brain cancer, with a 12% mortality rate, stands as the most lethal cancer among children. Out of the annual 4,000 diagnoses, 30% involve low-grade gliomas, with a troubling 10% progressing to lethal high-grade gliomas within 2 years. While surgery can often cure the remaining 90%, 50% necessitate adjuvant therapies like chemotherapy and molecular targeting, with radiation sparingly used. Alarmingly, 66% of survivors grapple with long-term health, mental, and cognitive issues due to treatment toxicity, emphasizing the urgent need for innovative biological therapies in pediatric glioma treatment.Our proposed strategy involves the acute analysis of freshly discarded tissue from brain cancer craniotomies, departing from the use of FFPE-preserved tissue, which may exhibit suboptimal quality. All patients provided informed consent, and by prospectively preserving viable tissue, we enable its utilization in ten or more subsequent analyses. These encompass cell-culture-based assays, viability assays, bulk, single cell, and single nucleus RNA sequencing, DNA sequencing, single-cell proteome analyses, flow cytometry, and immunofluorescence.Our findings highlight a direct correlation between surgical sample size and analysis feasibility. Bulk RNA sequencing proves valuable in acutely isolated surgical tissue and cryostored or flash-frozen tissue for identifying novel therapeutic targets in individual patients. Viable cells retrieved from cryopreserved tissue of the same patient allow ex vivo functional validation of therapeutic targets, eliminating the need for establishing patient-derived cell lines and xenografts, particularly challenging for low-grade gliomas. These findings establish the feasibility of a personalized bedside-to-bench approach, contingent on an adequate supply of surgical material.Expanding this approach to encompass all childhood brain cancer types, we plan to investigate using our comprehensive biorepository, housing viable tissue and single-cell suspensions from 122 childhood brain cancer patients. Among these, 51 cases yielded tissue for single-cell RNA-seq and viable cryoprotection. The most common tumor histology identified was low-grade glioma (32), primarily pilocytic astrocytoma (19), followed by embryonal (12) and ependymal tumors (4).In summary, we developed a standard operating procedure for analyzing fresh and cryopreserved surgical discards, identifying novel targets in childhood brain cancer tissue, and functionally validating them in patient-matched samples. This personalized medicine approach holds promise in circumventing the tedious development of patient-derived models. Citation Format: Anna Pronina, Emon Nasajpour, Caitlynn Tran, Troy Schouten, Eden Gibson, Jeffrey Nirschl, Laura Marie Prolo, Kelly Mahaney, Hannes Vogel, Claudia Katharina Petritsch. Biorepository yields potential for personalized medicine approach in pediatric brain cancer patients [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 5474.